X-ray spectra from convective photospheres of neutron 
stars 
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Abstract. We present first results of modeling convective 
photospheres of neutron stars. We show that in photo- 
spheres composed of the light elements convection arises 
only at relatively low effective temperatures ( < 3 — 5 x 10** 
K), whereas in the case of iron compositon it arises at 
^eff ^ 3 X 10^ K. Convection changes the depth depen- 
dence of the photosphere temperature and the shapes of 
the emergent spectra. Thus, it should be taken into ac- 
count for the proper interpretation of EUV/soft-X-ray ob- 
servations of the thermal radiation from neutron stars. 



1. Introduction 

Recent ROSAT observations of pulsars revealed that some 
of them emit thermal- like radiation in the soft X-ray range 
(Ogelman 1995). These observations stimulated further in- 
vestigations of neutron star (NS) photospheres responsible 
for the properties of the emitted radiation. 

One of important and hitherto untouched problems 
associated with modeling of the NS photospheres is the 
problem of convective energy transport which can affect 
the temperature distribution and the emergent spectra. 
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NS photospheres are much denser, p ^ 0.001 — 1 
g cm^"^, than those of usual stars. The increased densi- 
ties shift ionization equilibrium: the nonionizcd fraction 
grows with p at moderate densities, p < 0.01 — 0.1 g cm^'^, 
and sharply decreases at p > 0.1 — 1 g cm^'^ due to pres- 
sure ionization. As a result, zones of increased opacity 
(increased radiative gradient Vrad) and reduced adiabatic 
gradient Vad develop in the NS photospheres at temper- 
atures much higher than in photospheres of usual stars, 
which may cause convective instability at depths where 
the emergent spectrum is formed. So far the NS photo- 
spheres have been considered either for stars with very 
strong surface magnetic fields B ^ 10^^ — 10^^ G (see, 
for example, Pavlov et al. 1994) where existing convection 
theories are not applicable, or for 'nonmagnetic' (low-field) 
photospheres, B < 10^ — 10^ G, with high surface effective 



temperatures T^s > lO^ '^ - 10*^ K (Romani 1987), where 
the convection can hardly be expected. Since thermal NS 
radiation can be detected for Toff > 2 x 10^ K in the soft 
X-ray range, and for even lower temeratures, > 2 x 10* K, 
in the UV/optical range (Pavlov et al. 1995), the study 
of the convective NS photospheres is important for the 
proper interpretation of these observations. 

The convective flow in stellar photospheres is turbu- 
lent and imposes many complicated problems. A com- 
mon practice, which we follow here, is to use the phe- 
nomenological mixing-length theory with the traditional 
Schwarzschild criterion for convective instability, Vrad > 
Vad- This theory has been widely implemented and testi- 
fied. 

Details of our numerical calculations will be described 
elsewhere. Generally, we employ the complete lineariza- 
tion method for computing the photosphere models and 
include the convective energy transfer as described by Mi- 
halas (1978). 

2. Results and Discussions 

We calculated the model of nonmagnetic photospheres for 
different chemical compositions. Here we present examples 
for pure hydrogen, helium and iron compositions at the 
gravitational acceleration g = 2.43 x 10^"* cm s~^ which 
corresponds to standard NS mass M = 1.4Mq and radius 
i? = 10 km. The raditive opacities and equation of state 
for the iron composition were taken from the OPAL library 
(Iglesias et al. 1992). The results can be directly applicable 
to very old NSs with low magnetic fields (e. g., millisecond 
pulsars) . 

Our results show that, similar to the case of usual stel- 
lar photospheres, the convective energy transfer begins to 
play a role at lower surface temperatures when atoms are 
not fully ionized and the radiative opacities are strongly 
increased by contribution of the bound-free and bound- 
bound transitions. The increased opacities result in high 
values of Vrad — the depth dependence of the tempera- 
ture becomes steeper in order to transfer the energy flux 
throughout the photosphere. On the other hand, Vad in 
the dense partially ionized layers can be much smaller than 



its limiting value 0.4 for an ideal fully ionized or fully non- 
ionized gas (see, for example, Cox and Guili 1968). As 
a result, superficial convection zones in NS photospheres 
arise and, as in usual stars, they are associated with lay- 
ers of partial ionization. Our computations show that Vad 
can drop down to ~ 0.1 when the nonionized fractions are 
> 70%. 

The actual temperature gradient V in the convection 
zones satisfies the following relation: Vad < V < Vrad- 
Fig. 1 shows the temperature distributions in photo- 
spheres with and without allowance for convection. One 
can see that convection leads to more gradual profiles, in 
accordance with the above relation. Both in very surface 
layers, where Viad is too low, and in deep layers, where 
Vad is close to its maximum value, convection is absent, 
and the temperature profiles remain the same. 



Fig. 2. Spectra of outgoing radiation corresponding to the 
temperature profiles in Fig. 1. Dash-dotted curves are the 
blackbody spectra 7rB(refi). 



Fig. 1. Temperature profiles with and without allowance for 
convection (solid and dashed curves, respectively) for H, He 
and Fc photospheres {Tes = 1.8 x 10'', 4.1 x 10"* and 1.4 x lO'^ 
K, respectively). 

The convcctivc transfer affects not only the structure 
of photosphere but also the spectra of the NS thermal 
radiation (Fig. 2) because the temperature profiles are 
changed in the layers where the radiation escapes from. 
In particular, convection substantially (up to two orders 
of magnitide) lowers the flux from H and He photospheres 
at photon energies above the main photoionization edges, 
so that the high-energy spectral tails become softer. The 
spectra remain the same at low and very high energies 
since both shallow and very deep layers are not affected 
by convection. In the case of Fc composition, the convec- 
tive zone lies so deep that only a high-energy tail {E > 0.3 
keV) is affected. The effect of convection on the spectra 
disappears with increasing effective temperature (e. g., at 
Teff > 3 X 10" K for H, at Tes > 5 x 10"^ K for He and at 
Teff > 3 X 10^ K for Fe photospheres). 

The presented results correspond to the convective ef- 
ficiency l/H = 1. Acceptable values of this parameter can 
vary from 0.3 to 2.5. Bergeron et al. (1992) showed that 



higher efficiency enhances convection and smoothes the 
spectra of white dwarfs. The same effect can be expected 
in the NS photospheres and, consequently, the convection 
there may develop at higher effective temperatures than 
in the models presented. 

Convection in NS photospheres is important because it 
can mix the material in convective zones, bringing heav- 
ier elements from bottom to surface layers which would 
otherwise contain only light elements due to gravitational 
stratification. Our calculations show that this may happen 
in cold NSs with T^s < (3 - 5) x 10* K, whereas radiating 
layers of hotter NSs can be expected to consist mainly 
of hydrogen and helium. Since the presence of convec- 
tion softens the high-energy tails of the spectra, this effect 
should be taken into account for the proper interpretation 
of EUV/soft-X-ray observations of thermal radiation from 
NSs. 
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